
qui{

use "$datapath\001_nlsy_child_mother_cog.dta", clear

replace wage_hour_beforebirth_mom = . if wage_hour_beforebirth_mom==0

** Delete top bunching
drop if hours == 2080

 
histogram hours if hours>0 & hours<=100, width(2) frequency  ///
       scheme(s1mono) legend(off) ///
	   ytitle("Number of Mothers in the Sample") ///
	   xtitle("Average hours working during three first years of child") ///
	   fcolor(gs9) lcolor(black)
graph export "$figpath/002_BunchingPDF_hours.pdf", as(pdf) name("Graph") replace	   

}